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Abstract 

^, 
CL|| We analyze the light cone (loffe) time structure of the gluon dis- 

o ' tribution function in the double logarithmic approximation. We show 

O . that due to QCD evolution loffe equation is modified. The characteris- 

tic light cone time of the gluons does not increase as fast with increas- 
ing energy (decreasing Bjorken x) as predicted by the parton distribu- 
^ ■ tions exhibiting Bjorken scaling due to the increase of the transverse 

ixLomenta of the gluons in the DGLAP ladder. 



1 Introduction 

It was first observed in Quantum Electrodynamics |^ that photon emission in 
the scattering process of electrons propagating through a medium occurs over 
distances increasing with energies. It their seminal paper Gribov, loffe 
and Pomeranchuk have demonstrated that at high energies large longitudinal 
distances, which are now usually referred to as coherent length distances, Ic, 
become important in the virtual photon-nucleon scattering when considered 
in the rest frame of the target. (One can also formulate this statement in 
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the frame independent language by using light-cone time along the reaction 
axis.) In the following years loffe [Q built an explicit connection between the 
deep inelastic scattering amplitude and the space-time representation of the 
correlator of the electromagnetic currents. He analyzed the first deep inelastic 
data from SLAC assuming the exact Bjorken scaling and demonstrated that 
the longitudinal distances in the the Bjorken limit are growing as 

h ex ^, (1) 

niN X 

where tun is the target nucleon's mass. A natural implication of the increase 
of longitudinal distances was the parallel observation by Gribov 0] that large 
shadowing effects should be present in the deep inelastic scattering off nuclei 
at small enough x when /^ ^ 2_R^, leading to the saturation behavior of 
F2A0cA^/^QHn{l/x). 

Qualitatively one can understand the pattern of Eq. (0) by considering 
the energy denominator for the transition of the virtual photon to some 
excited Fock state consisting of quarks and gluons with the invariant mass 
M^. Indeed the uncertainty principle tells us that the lifetime of this state 
in the nucleon rest frame is the inverse of the energy non-conservation in the 
7* -^ "M" transition and hence ~ 2qo/{M'^ + Q'^), which for M^ comparable 
to Q^ leads to Eq. (^. Bjorken P has build an explicit picture of the space 
time evolution of the quark-antiquark pair in the aligned jet model of DIS 
for the transition 7* — >■ qq and demonstrated that even though the pair is 
produced at zero transverse separation it reaches hadron size scale over the 
time comparable to the coherent length /^ of Eq. (0). 

Further studies of Ic were performed in in the aligned jet model. It 
was demonstrated that in this model the average value of M^ is (M^) ^ Q^, 
leading to Ic = ^m x ' Similar numbers were obtained in ^ by explicit 
consideration of the loffe representation assuming that the structure function 
exhibits Bjorken scaling and that the nucleon's structure function behaves 
as F2n{x) —>■ const in the x — >■ limit. 

More recently in a number of papers the coordinate representation of 
nucleonic and nuclear correlators of electromagnetic currents was considered 
using the current parameterizations of parton densities and employing the 
Bjorken scaling approximation for the loffe expressions 0, §]. The essential 
distances in the coherent diffractive interactions with nuclei which lead to the 
leading twist shadowing were investigated in B and coherence distances were 



found to be of the order of ~ l/2mNX for moderate Q^ but were decreasing 
with increasing Q^. 

The aim of this paper is to investigate how the scahng violation which 
is especially strong in the region of small x affects previous conclusions for 
the value of Ic- Qualitatively we can expect that due to production of inter- 
mediate states with a large invariant mass in the ladder-type kinematics the 
coherence time Ic should be substantially reduced. As a first step we will use 
a double logarithmic approximation which emphasizes the ladder kinematics. 
In Sect. 2 we will consider the double logarithmic evolution with fixed strong 
coupling constant and in Sect. 3 we will analyze the case of running coupling 
constant in the evolution equations. 

In Sect. 4 we will study how the new value for Ic generated by the double 
logarithmic evolution affects the onset of non-linear effects such as saturation 
[^, ^ |rT|, |T^] in a nucleus. To do that we will have to estimate at what 
energies the coherence length becomes larger than the nuclear diameter 2i?^. 
As we will observe the double logarithmic evolution slows down the onset of 
these phenomena with energy. 

We want to analyze the structure of the gluon distribution in a proton 
in terms of the loffe time variable u = p ■ x, where p^ is the target nucleon's 
momentum and x^ is the parton's coordinate. To avoid confusion throughout 
the paper we will write z for the Bjorken x variable. The gluon distribution 
function can be represented as a Fourier transform in light cone time u as 
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with u = p ■ q/rriM and Q^ = —q^ the photon's virtuality. Here the function 
f{u,v) in general also depends on the transverse momentum squared q^. 
Inverting the Fourier transform of Eq. (^ we write 

f{u,h') = / ■-— e'''"' zG{z,2mNJyz) — complex conjugate. (3) 

JO 27r 

Therefore, to determine the coherence length l^, or, equivalently, the char- 
acteristic loffe time u of the gluon distribution we have to find the function 
f{u, v) using the inverse Fourier transform of Eq. (^ and extract the typical 
values of u out of it. 



2 lofFe Time Evolution of the Double Log- 
arithmic Structure Function: Fixed Cou- 
pling Case 

We begin by considering the gluon distribution function of a hadron in the 



double logarithmic approximation. The DGLAP equation [|14| for the gluon 



structure function in the double logarithmic approximation is 



where the gluon-gluon splitting function at small z is given by 

iGciz) ^ -^. (5) 

To solve Eq. (^ we for simplicity choose the following initial condition at 
some not very large value of photon's virtuality Ql 

Giz,Ql) =6iz-zo). (6) 

In this section we will consider the "toy model" case when the strong coupling 
constant is fixed. Then the solution of Eq. (|^) is given by 

cisNc 1 
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zG{zM')= ^A-] ^ , (7) 
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where the integral over n runs along a straight line parallel to imaginary axis 
to the right of all the singularities of the integrand. Using Eq. (|^) in Eq. (|^) 
we obtain 
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— complex conjugate. (8) 

Performing the z integration in Eq. (H) yields 

n n — 1, 
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where the integration over n can be approximated by the saddle point method 
around 



to give 
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Eq. (|Tl]) gives the Fourier image in loffe time of the gluon structure function 
in a proton. Our goal is to determine which values of the variable u dominate 
in Eq. (§) with the function /(u, z/) given by Eq. ([TT|) . Most of the z and 
u dependence in Eq. (pA| ) is driven by the exponential function in it. To 
determine the characteristic values of u we have to find which u gives the 
maximum of the function in the exponent. One can readily see that the 
maximum is reached at 

« c = - \-- (12) 

Another way of determining the typical value of u is by finding the median 
of the integral over u in Eq. (|]) 0. Our numerical estimates showed that 
the discrepancy between that method and the one presented above is not 
significant in the double logarithmic approximation and thus Eq. (|12D does 
give a good estimate of the typical u in the gluon distribution function. 

One can see that the gluon smearing of Eq. ( ]T^ ) is not as large as one 
would naively expect from the loffe equation where in the Bjorken scaling 
limit there assumed to be no z/ dependence in the integrand of Eq. (|^) . There 
the spread of partons is limited only by the Fourier exponent in Eq. (|^) 
yielding the maximum possible value of u to be 

''^max — ~ (13) 



which translates itself into the well-known coherence length of the gluon in 
the proton's rest frame Icoh = l/2mNZ. Recalling that u = Q'^/2mNZ one 
can see that the time in Eq. (^) grows much slower with decreasing z than 
the time given by Eq. (p!3D. 

To understand what happens it is instructive to rewrite Eq. (|T^ ) in terms 
of the x_ coordinate and g_ momentum of the last gluon in the ladder. 
Since we are working in the infinite momentum frame where p^ = (p+,0,0) 
we write u = p ■ x = p+x_ and v = p ■ q/rrijq = p+q_/mN which leads us 
to the following expression for the light cone lifetime of the last gluon in the 
ladder 

.;_ . ^J^ . 41^, (14) 



/Q^Zo2p+q^ \IQIzqS 
where we introduced the center of mass energy s = {p + q)"^ ~ 2p^q_. On the 



other hand in a frame where several evolution rungs start in the probe [12 



we can write down the light cone lifetime of a (longitudinally) soft gluon as 

x*_ = %. (15) 

To estimate the typical g^ of a gluon in the DLA DGLAP evolution taken at 
a given fixed value of g_ (or, equivalently, s) we first note that in the saddle 
point approximation Eq. (^ leads to 



zG{z, Q') oc exp hJ^ In S In ^ J • (16) 

Substituting Q"^ ^ g^ and z ^ q^ /2p^q_ we obtain the typical transverse 
momentum squared of the gluons 



t)j, - ^Qlzo2p+q^ ^ \/QlzoS, (17) 

which, after being substituted into Eq. (|15]) gives us Eq. ([T^). Now the 
difference between loffe times in the distribution function exhibiting Bjorken 
scaling and the double logarithmic gluon distribution becomes apparent. As 
we increase the center of mass energy s by increasing q_ the typical transverse 
momentum of the distribution function without any QCD evolution would 
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remain constant, while in the case of double logarithmic evolution it would 
increase with energy (see Eq. i\rf\j ) leading to shorter light cone lifetimes 
in this latter case. The diffusion of the typical transverse momentum with 
energy is so fast in our Eq. (|l^ because in this section we considered the toy 
model of fixed coupling DLA DGLAP equation. 

3 lofFe Time Evolution of the Double Loga- 
rithmic Structure Function: Running Cou- 
pling Case 

Let us repeat the calculation of the previous chapter for the running coupling 
case. The running coupling constant in Eq. (^ will be taken below at the 
one-loop level 

with the beta function b = — ^ — -. 

The solution of the double logarithmic DGLAP equation for gluon dis- 
tribution function (^ with the splitting function of Eq. (^) and the initial 
conditions given by Eq. (H) is 



dn fzoy-^ /ln(QVA2 



Nc 
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where, just like in Eq. (|^) the integral over n runs along a straight line 
parallel to imaginary axis to the right of all the singularities of the integrand. 
Substituting Eq. ( [I9| ) into Eq. (|^) we write for the loffe time representation 
of the gluon structure function ^ 



., , , dz _-^ f dn /-2o\" ^ / \n(2iymNz/A 
f{u,iy) = / — e ™" ' ' ^ ' ' 
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— complex conjugate. (20) 



Integration over z in Eq. pO| ) is rather complicated but can be simplified if 
with the leading logarithmic accuracy we substitute z by Iju in the ratio of 
the logarithms in it. The rest of the integral can be done easily yielding 
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The expression in Eq. ( |2TD can be evaluated by the saddle point method. 
The position of the saddle point is given by 



no 



iVc ( ln(2z/miv/MA2) 
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and Eq. ( pT]) becomes 
/(m,z^) = 
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Eq. ( p3D yields the distribution in light cone loffe time of the gluon structure 
function in a proton. Similarly to the previous section in order to find the 
characteristic value of u we need to find the maximum of the function in 
the power of the exponent in Eq. (^). Defining L = ln(2z/r?7,7v^o/^^) and 
Lq = ln(2i/o'TiArZo/A^) together with ( = \ia{uzo)/L allows us to rewrite the 
function under the square root in the exponent of Eq. (ESF) as 



/i(C) = C^ln 



Ln 



(24) 



For L ^ Lq the maximum of Eq. ( |24D is reached at the value of C that can 
be very well approximated by the formula 

1 



c* 



1 



InL/Lo 



(25) 



which is found from the usual extremum condition 

h'{C) = In — + ln(l - C*) - -^ ^ In — - ^ -^ = 0. (26) 



Eq. (|25D can be rewritten in terms of the loffe time and u as 

l/jn(2i;™jvfo/^A 

^ ~ - ^2 (27) 



zo V A2 
or in terms of x_ and g_ as 



2g_ (2p+q^ZQ 



A2 V A2 



V in(Qj5/A2) y 



A2 Ia2 



(28) 



To understand this result we, similarly to the previous section first note 
that in the saddle point approximation the gluon distribution of Eq. (^) is 
proportional to 



.G,.,e^)ocexp|2^^,„(|^)l„^| (29) 

with the typical transverse momentum for fixed g_ given by 
/ 2\ ^ ^2 /2p+g_^\ I in(Q^/A2) J ^ ^2 /"^y" Imq^T^J _^3q^ 



A2 

Substituting Eq. ( |30D into Eq. (p!5| ) readily yields us Eq. (pSl). Thus again 
the DGLAP evolution makes transverse momenta in the gluon distribution 
diffuse towards larger values slowing down the spreading of the gluons in the 
longitudinal x_ direction as compared to the standard estimate of Eq. ([T3|) . 
The growth of the transverse momenta with s described by Eq. (RUl) is slower 



than any positive power of s. Thus it is slower than (g^y ~ ^/s of the fixed 



coupling case (see Eq. (|17D). Slower growth of the transverse momentum in 
the running coupling case leads to longer light cone coherence times for the 
gluons than in the fixed coupling case, as follows from Eq. (|T5|). Nevertheless 
the coherence time is still shorter than the maximum limit of Eq. (|13|). 
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4 Discussion 

To illustrate our results we are going to plot the gluon coherence lengths 
produced by different evolution scenarios discussed above in the rest frame 
of the proton as functions of center of mass energy s. By doing so we would 
also address the question posed in the Introduction: at which values of z does 
the DGLAP evolved gluon distribution of a nucleon in a nucleus reach the 
lengths comparable to the nuclear diameter, so that the non-linear effects 
involving multiple rescatterings and mergers of gluon ladders would begin 
taking place? 

In order to plot the coherence length of the gluons in the rest frame of 
the target proton or nucleus we use the simple relationship u = l/2mj^lcoh 
and employ Eqs. (0) and (28) for the fixed and running coupling cases 



correspondingly. We use A = 0.3 GeV and 2;o = 0.1. In order for the diffusion 
of transverse momenta to start at the same initial value we have to choose 
different values of the initial virtuality Qo for the fixed and running coupling 
cases: we use Qo = 4 GeV in Eq. (0) and Qo = 0.8 GeV in Eq. (|2|). By 
doing so we make the initial transverse momentum squared of both evolutions 
equal to q^^ ^ 20 GeV^ at s = 200 GeV^. 

The plot of Icoh as a function of the center of mass energy s is depicted 
in Fig. ||. The medium-thick line corresponds to the fixed coupling DLA 
DGLAP evolution [IJ^) while the thick line corresponds to the running cou- 
pling case {l*c)- The top solid line depicts the upper limit on the coherence 
length given by Icoh = l/2mNZ, where z = qfnit/^ ^o insure that this length 
is really maximum possible. Finally the horizontal line corresponds to the 
diameter of some sample large nucleus taken here to be equal to 2R = 14 fm. 
One can see that due to DGLAP evolution the coherence lengths of gluons 
do not grow as fast with energy s as predicted by the well-known estimate of 
l/2?Tijv-2 (upper line). For the case of running coupling the coherence length 
is roughly a half of the maximum throughout the whole region of s consid- 
ered. Appropriately, in the case of a nuclear target, the coherence length of 
the partons in each individual nucleon reaches the nuclear diameter at higher 
values of s than expected from Eq. (plSf). In the fixed coupling case the nu- 
clear diameter is reached at the energies beyond those shown in Fig. [1|, at 
about s ~ 30000 GeV^. In a more realistic case of running coupling DGLAP 
evolution the crossover happens slightly above the one predicted by Eq. (p!3|), 
leading to about 30% higher s required for the gluon's coherence length to 
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Figure 1: Distances characterizing the spread of gluons in the rest frame of 
the proton. Horizontal hne denotes the diameter of some large nucleus which 
is taken here to be 14 fm. The top solid line corresponds to the upper limit of 
the Eq. (|T3p, the lower (medium-thick) line corresponds to the fixed coupling 
DLA DGLAP evolution of Eq. (0) and the middle (thick) line depicts the 
running coupling evolution of Eq. (^). 



become comparable to the nuclear radius (see Fig. [l|). To summarize, based 
on the example of the double logarithmic approximation we conclude that 
DGLAP evolution towards higher Q^ slows down the onset of non-linear 
effects in the nuclear wave function. 

From the picture presented in Fig. |I] one concludes that much smaller 
than naively expected values of z (for the same large Q^) are needed in 
the nuclei for the parton distributions of the individual nucleons to start 
overlapping making non-linear effect related to nuclear shadowing and parton 
saturation possible. However one must keep in mind that the spread of the 
gluon distribution in Fig. [^ is generated by the DLA DGLAP evolution, that 
is along with going towards smaller values of z the evolution also moves 
towards larger values of transverse momentum g^ as could be seen from Eqs. 
(O) and (pOl). Going towards smaller z tends to increase light cone times 
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and corresponding coherence lengths in the target's rest frame. At the same 
time increasing q^ tends to decrease hght cone times, moving the evolution 
away from the saturation region of moderate q^ ~ Ql. Thus even though 
DGLAP evolution seems to slow down the growth of the light cone times 
as functions of z as could be seen from Eqs. (|14D and (0), this is entirely 
due to the fact that DLA DGLAP evolution moves the distributions toward 
larger q^ pushing it away from the saturation region and non-linear effects. 

Of course the gluon distribution generated by the double logarithmic 
evolution equation describes the data well only in a rather narrow kinematic 
region of small x and large Q^ . Therefore our analysis of light cone loffe time 
structure of the DLA gluon distribution is also limited to this kinematic 
region. Further studies involving more realistic parton densities beyond the 
double log approximation are necessary in order to fully quantify the effect 
of scaling violations on the light cone structure of distribution functions. 

Acknowledgments 

The authors are very much indebted to Lonya Frankfurt for several insightful 
discussions. M.S. thanks the Department of Energy's Institute for Nuclear 
Theory at the University of Washington for its hospitality and the Depart- 
ment of Energy for partial support during the time this study started. The 
work of Yu.K. was supported in part by the U.S. Department of Energy under 
Grant No. DE-FG03-97ER41014. The research of M.S. was sponsored in part 
by the U.S. Department of Energy under Grant No. DE-FG02-93ER40771. 

References 

[1] L. Landau, I. Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 92, 535 (1953). 

[2] V.N. Gribov, B.L. loffe, I.Ya. Pomeranchuk, Sov. J. Nucl. Phys. 2, 549 
(1966); Yad. Fiz. 2, 768 (1965). 

[3] B. L. loffe, JETP Lett. 9, 163 (1969); 10, 143 (1969); Phys. Lett. 30B, 
123 (1969). 

[4] V. N. Gribov, Sov. Phys. JETP 30, 709 (1970), Zh. Eksp. Teor. Fiz. 57, 
1306 (1969). 

12 



[5] J. D. Bjorken, Electroproduction at Very Small Values of the Scal- 
ing Variable, contribution to the Festschrift of Jack Steinberger, K. 
Kleinknecht and T.D. Lee (editors). 

[6] L. Frankfurt and M. Strikman, Phys. Kept. 160, 235 (1988). 

[7] G. Filler, M. Vanttinen, L. Mankiewicz, W. Weise, K. Eskola, Nucl. 
Fhys. A663, 328 (2000). 

[8] S. Hofmann, L. Mankiewicz, A. Schafer, Z. Fhys. A359, 157 (1997); V. 
Braun, F. Gornicki, L. Mankiewicz, Fhys. Rev. D51, 6036 (1995). 

[9] L. Frankfurt, V. Guzey and M. Strikman, J. Fhys. G 27, R23 (2001). 

[10] L. V. Gribov, E. M. Levin, and M. G. Ryskin, Fhys. Rep. 100, 1 (1983); 
A.H. Mueller, J.-W. Qiu, Nucl. Fhys. B268, 427 (1986). 

[11] L. McLerran and R. Venugopalan, Fhys. Rev. D 49, 2233 (1994); 49, 
3352 (1994); 50, 2225 (1994); Yu.V. Kovchegov, Fhys. Rev. D 54, 5463 
(1996); 55, 5445 (1997); J. Jahhan-Marian, A. Kovner, L. McLerran, 
and H. Weigert, Fhys. Rev. D 55,5414 (1997). 

[12] I. I. Balitsky, Nucl. Fhys. B463, 99 (1996); Yu. V. Kovchegov, Fhys. 
Rev. D 60, 034008 (1999). 

[13] B. L. loffe, V. A. Khoze, L. N. Lipatov, Hard Processes, vol. 1, Elsevier 
Science Fublishers, (1984). 

[14] Yu. L. Dokshitzer, JETF 73, 1216 (1977); V. N. Gribov and L. N. 
Lipatov, Sov. J. Nucl. Fhys. 15, 78 (1972); G. Aharelli and G. Farisi, 
Nucl. Fhys. B126, 298 (1997). 



13 



